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This is the second in a series of papers describing the ALPAK system 
for nonnumerical algebra on a digital computer. The first paper, Ref. 1 , is 
concerned with polynomials in several variables and truncated power 
series with polynomial coefficients. This paper is concerned with rational 
functions of several variables and truncated power series with rational- 
function coefficients. A third paper, Ref. 3, will discuss systems of linear 
equations with rational-function coefficients. 

The ALPAK system has been programmed within the BE-SYS-4 
monitor system on the IBM 7090 computer, but the language and concepts 
are machine-independent. Several practical applications arc described in 
Ref. 1. 

This paper is divided into five sections. The fi,7st deals with basic con- 
cepts, the second defines canonical forms, and the third describes ALPAK's 
greatest common divisor algorithm. These three sections do not presuppose 
any knowledge of computers or computer programming. Section IV de- 
scribes the use and the implementation of the algebraic operations relating 
to rational functions of several variables and truncated power series with 
rational-function coefficients. The reader of this section is assumed to be 
familiar with the basic concepts of computer programming and with Ref. 1. 
Finally, Section V discusses very briefly some of our plans and hopes for 
the future. 
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I. SUMMARY OF THE AVAILABLE RATIONAL-FUNCTION OPERATIONS 

l.i Introduction 

The ALPAK system is a programming system for performing routine 
manipulations of algebraic expressions on a digital computer. The 
system operates on rational functions of several variables and on trun- 
cated power series in several variables with rational functions of several 
other variables as coefficients. It is capable of performing the operations 
of addition, subtraction, multiplication, division, substitution and 
differentiation. In the present version of the program the coefficients 
of the rational functions are integers, but the change to coefficients 
from any other integral domain can be made without major program 
reorganization. ALPAK is also capable of solving systems of equations 
linear in certain variables with coefficients which are rational functions 
of other variables (see Ref. 3). The ALPAK system as described in this 
paper has been programmed for the IBM 7090 computer. 

1.2 Input-Output 

Rational functions can be entered into the machine from punched 
cards, and the output can be printed and /or punched. The polynomial 

P(x,y,z) = 8xy 2 z + 2xifz 2 - 10x*yz z 

can be entered into the machine by punching the following array of 
coefficients and exponents one term per card: 

8 1,2,1 

2 1,2,2 
-10 3,1,3 
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The zero coefficient is an end-of -polynomial signal. The subroutine 
which reads these cards must have access to a polynomial format state- 
ment (previously read from cards) containing the names of the variables 
and the number of bits to be allocated for the exponents of each. One 
full word (35 bits plus sign) is allocated for the coefficient of each term, 
thus permitting coefficients up to 2 35 — 1 in magnitude. Rational func- 
tions which are not polynomials are entered by punching the numerator 
and denominator polynomials successively and calling the rational- 
function reading subroutine. 

The punched and printed output consists of arrays of coefficients and 
exponents similar to those that are accepted as input. Input and output 
are accomplished by simple commands such as 



read format statement FMT 
read rational function R 
print rational function R. 



RFNRDF 


FMT 


RFNRDD 


R,FMT 


RFNPRT 


R 



1.3 An Example of the ALPAK Language 

The simplicity of ALPAK programming is illustrated by the following 
example. Suppose rational functions A, B, C, and D and a format state- 
ment FMT have been punched on cards, and we wish to compute and 
print the rational function 

F = (A*B/C) + D, 

where the asterisk, denotes multiplication. 
The required program is 

begin (reserve 10000 words of storage 
for data and working space) 

read polynomial format statement FMT 
from cards 

read polynomial A from cards 

read polynomial B from cards 

read polynomial C from cards 

read polynomial D from cards 

replace F by A *B 

replace F by F/C (C must not be 
zero) 

replace F by F + D 

print F 

go to END JOB. 



RFNBEG 


10000 


RFNRDF 


FMT 


RFNRDD 


A,FMT 


RFNRDD 


B,FMT 


RFNRDD 


C,FMT 


RFNRDD 


D,FMT 


RFNMPY 


F,A,B 


RFNDIV 


F,F,C 


RFNADD 


F,F,D 


RFNPRT 


F 


TRA 


ENDJOB 
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II. CANONICAL FORMS 

2.1 Introduction 

All rational functions stored by the program are kept in a unique 
canonical form which is the subject of this section. The read routines 
place the input functions in canonical form and all operations leave 
their results in canonical form. The uniqueness of the canonical form 
ensures that two equal rational functions with the same format are 
precisely identical in storage and that, in particular, zero is uniquely 
represented. 

2.2 Polynomial Canonical Form 

A polynomial is always represented as an ordered list of its nonzero 
terms. It is convenient to order the terms according to the magnitude 
of the first exponent, and to order those terms having the same first 
exponent according to the magnitude of the second, etc. The order of 
the variables is the order in which they appear in the format statement. 

2.3 Rational-Function Canonical Form 

A rational function is represented as an ordered pair of polynomials, 
namely its numerator and denominator respectively. These must be in 
polynomial canonical form, and they must be relatively prime. In addi- 
tion, the sign of the numerator must be chosen so that the first term of 
the denominator is positive. 

III. THE GREATEST COMMON DIVISOR ALGORITHM 

3.1 Introduction 

Since rational functions in canonical form must have numerator and 
denominator relatively prime, the ALPAK program must be capable 
of finding the greatest common divisor (G.C.D.) of polynomials in 
several variables. This is the essential ingredient in the extension of 
ALPAK from polynomials to rational functions. Since each rational- 
function operation must leave its result in canonical form, the G.C.D. 
operation will be performed very frequently in most programs involving 
rational functions. 

Let «i , a 2 , • • • , a n be a set of nonzero polynomials. A G.C.D. of 
fli , <h , • • • , a» is defined to be a polynomial g such that 
(i) g divides each of a,\ , flfc , • ■ • , a n ; and 
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(ii) any polynomial g' that divides each of Oi , 02 , • • • , o B also divides g. 
We denote a G.C.D. of a, , a 2 , • • ■ , a n by (a x , a 2 , • • • , o») . Since every 
polynomial has a decomposition into primes that is unique up to sign, this 
definition implies that a G.C.D. is unique up to sign. In the special case 
of integers, the positive value is often referred to as the G.C.D. 

The next three subsections discuss the Euclidean algorithm for in- 
tegers, ALPAK 's generalization of it for polynomials, and some special 
strategies which make the latter more effective. The final subsection 
attempts to present a balanced picture of the present capabilities of the 
ALPAK algorithm. 

Algebraic background relevant to the following discussion can be 
found in Chapter 1 of Ref. 2, or almost any other algebra text that 
treats polynomial rings. 

3.2 The Euclidean Algorithm 

The G.C.D. of a set of n nonzero integers can be obtained by a series 
of pairwise computations, because 

(fti , • • ■ , a n ) = (((•• • ((01 ,a 2 ),a 3 ), • • • ),On_i),o»). (1) 

The Euclidean algorithm obtains the G.C.D. of two nonzero integers a 
and b. Without loss of generality we can assume that both are positive 
and that a ^ b. By the division algorithm we can write 

a = qb + c (2) 

with 

^ c < b. (3) 

If c = 0, then b divides a, so (a,b) = b. Otherwise the common divisors 
of a and b are the same as those of b and c, so (0,6) = (b,c). Since 
/ ; _|_ c < a + b, the process terminates in a finite number of steps. 

3.3 The ALPAK G.C.D. Algorithm 

We shall consider a polynomial in v variables as a polynomial in one 
variable, to be called x, with coefficients from the integral domain of 
polynomials in the remaining v - 1 variables. We shall represent these 
v - 1 variables by the vector y. If p(x,y) is such a polynomial, then 
d x (p) denotes the degree in x of p. 

Now let ft and b be a pair of nonzero polynomials. We shall present 
an inductive definition of the ALPAK algorithm, to be called POLGCD, 
for obtaining their G.C.D. Let v' be the number of variables in the pair. 
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If v' = 0, then a and b are both integers and the Euclidean algorithm 
is used. Assume POLGCD works for v' < v. We shall define it for v' = v. 
To begin, we write 

a(x,y) = a T (y)x T + a r . 1 (y)x r ~ 1 + ■ ■ • + a (y) 

(4) 
b(x,y) = b„(y)x 8 + b^(y)ar l + • • • + My). 

Our first task is to rewrite this as 

a(x,y) = xj(y)a'(x,y) 

8 (**' 

b(x,y) = x g(y)b'(x,y) 

where a' and b' are primitive in x; that is, neither is divisible by x or by 
any polynomial independent of x except =1=1. Clearly, a and /3 are the 
largest integers such that x* divides a(x,y), and x divides b(x,y); while 
f(y) is the G.C.D. of the nonzero o,(y), and g(y) is the G.C.D. of the 
nonzero bj(y). Since the a, and the by depend on fewer than v variables, 
our induction hypothesis implies that POLGCD will obtain / and g. 
Next we observe that 

(a,b) = (x a ,x )(f,g)(a f ,b'). (G) 

(The proof of this depends on the fact that a' and b r are primitive.) It 
is obvious that 

(*V) = * 7 (7) 

where 7 is the smaller of a and /3. Since / and g depend on fewer than v 
variables, our induction hypothesis implies that we can use POLGCD 
to obtain (f,g). 

We shall now define a subalgorithm, to be called PRMGCD, for 
obtaining the G.C.D. of the primitive polynomials a'(x,y) and b'(x,y). 
To begin, we write 

a'(x,y) = a m '(y)x m + a m ^'(y)x m - 1 + • • • + a Q '{y) 

(o) 
b'(x,y) = b n '(y)x" + b„- l '(y)x n - 1 + ••• + W(y), 

where aj ,a Q ' ,b„' ,b ' are all nonzero. Without loss of generality we can 
assume that d x (a') ^ d x (b'). If d x (b') = 0, then V = b ' = ±1,* so 
(a',b') = 1. Otherwise, we use POLGCD to compute 

h = (a m ',b n '). (9) 

Then we form 

* Here we have used the definition of primitivity which is giveu following (5). 
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iu which the bracketed fractions are polynomials. If c = 0, then a' = 
±6';* so (a',b') = 6'. Otherwise, let c'(x,y) be the primitive part of 
c(x,y), defined as in (6). Then the common divisors of a' and b' are 
the same as those of b' and c'* so {a',b') = (b',c'). (Note the similarity 
of this algorithm to the Euclidean algorithm, which was defined in the 
preceding subsection.) By construction 

d.(c') ^ d,(c) < d.(a'). (11) 

Hence 

d.(b') + d x (C) < d x (a') + d x (b'), (12) 

and so the process terminates in a finite number of steps. 

3.4 Special Strategies 

From a practical point of view the POLGCD algorithm leaves much 
to be desired. As the degree in x is reduced, the coefficients grow. This 
phenomenon is vividly illustrated by the following example. It is desired 
to find the G.C.D. of the primitive polynomials 

-39a: 4 + 125a: 3 - 15/ - 135a; - 44 , 

-12./ - 89/ + 192./ - Qx - 85. { } 

The successive pairs of primitive polynomials produced by the 
PRMGCD algorithm arc: 

-12x 4 - 89a; 8 + 192/ - Gx - 85 
- 1657a: 3 + 2556a; 1 + 462a; - 929 

178145/ - 312600a; 2 - 1206a; + 140845 
-1657a; 3 + 2556/ + 462a; - 929 

[-1657/ + 2550/ + 462x - 929 
5219965/ - 6692054a: - 5656955 

2253497062/ - 6961950605a; - 4849347485 
5219965/ - 6692054a; - 5656955 



To proceed farther would require double-precision coefficients, which 
are not now available in ALPAK. For polynomials in many variables 
this problem is even more acute. If POLGCD were programmed to 
handle coefficients and exponents of arbitrary size, the time it would 

* Here we have used the definition of primitivity which is given following (5). 



(14) 
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require to find the G.C.D. of two general polynomials in n variables, 
each having degrees d\ , • • ■ , d n respectively and each having coeffi- 
cients of modest size, can be shown to be proportional to 

, (c dl ) 

( 15 ) 

where c tt 2.6 is the square of the Fibonacci ratio, §(1 + V5). 

It is apparent that any real solution to the problem of coefficient 
growth would require a fundamentally different algorithm. However, 
many of the G.C.D. problems which arise in practice exhibit special 
properties which can be exploited. 

The most important of these special properties is variable independ- 
ence. If one of the inputs to POLGCD is independent of one or more of 
the variables, then the other can immediately be broken into subpoly- 
nomials, and we obtain a set of subproblems each involving only those 
variables which both of the original inputs depend on. For example, it 
is clear by inspection that the G.C.D. of the pair 

22 (2a; 4 - 17a; 3 + 65a; 2 - 144a; + 72) 
- 3y\2x 3 - 17a; 2 + 66a; - 72) 

6a- 4 - 41a: 3 + 104a; 2 - 116a; + 48 

is equal to the G.C.D. of the triple 

2a; 4 - 17a; 3 4- 65a; 2 - 144a; 4- 72 

2a; 3 - 17a; 2 + 66.r - 72 

6a; 4 - 41a; 3 + 104a; 2 - 116a: 4- 48, 

and the POLGCD algorithm will discover this provided that x is the 
last variable in the format statement. Unfortunately POLGCD does 
not now optimize the order of the variables. If the preceding example is 
attempted using the variable ordering (x,y,z), disaster ensues. The in- 
puts are viewed as 

4za; 4 - (6?/ 2 4- 34z)a; 3 + (51?/ + 130z)a; 2 
- (198?/ 4- 288z)a; + (216?/ + 144z) 

6a; 4 - 41a; 3 + 104a; 2 - 116a; 4- 48, 

and it is easily seen that both are primitive in x. The next two pairs of 
primitive polynomials produced by PRMGCD are 
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6.r 4 - 41.x- 3 + 104.x- 2 - 116a; + 48 

(18/ + 20z).r 3 - (153/ 4- 182z)a- 2 

+ (59/ + 4522) x - (648/ + 3362 ) 

(45/ + 682). x- 3 - (423/ -\- 158z).r. 2 

+ (450/ - 76z).x + (216/ + 240*) 

(18/ + 20z).x 3 - (153/ 4- 182z).x 2 

4- (59/ 4- 452z)z - (648/ -\- 336s). 



The variable independence has now been lost, and the subsequent 
pairs will have progressively higher degrees in y and z, and progressively 
larger coefficients. 

Both POLGCD and PRMGCD test their inputs to see whether 
either divides the other. Since a given G.C.D. problem may involve 
many recursive calls to POLGCD and PRMGCD, this strategy pays 
frequent dividends. POLGCD and PRMGCD also make full use of 
the fact that the G.C.D. of the set of terms of a polynomial, and simi- 
larly the G.C.D. of a monomial and a polynomial, can be computed 
simply and directly. 

Finally, we remark that the PRMGCD process is terminated as soon 
as the degree in x of either input is zero or one. A primitive polynomial 
of degree zero is obviously equal to ±1, while a primitive polynomial 
of degree one is irreducible. At the last variable level the PRMGCD 
process is terminated as soon as the degree in x of either input is three 
or less. A quadratic polynomial can be factored, if it is reducible, with 
the aid of the quadratic formula. A reducible cubic must have at least 
one rational root. A simple change of variable produces a related cubic 
which must have at least one integral root, and it is easy to test for this 
numerically. 

3.5 Concluding Remarks 

As we have already stated, the G.C.D. operation is the essential 
ingredient in the extension of ALPAK from polynomials to rational 
functions. The weakness of the ALPAK G.C.D. algorithm is apparent 
from (15). Its strength lies in the fact that most G.C.D. computations 
which arise in problems of practical interest have a degree of immunity 
from that formula because of their special structure. 

As an example we wish to mention the problem of a single-server 
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queue with feedback. The computation of the first two moments of the 
total time is outlined in Section II of Ref. 1 and in the Appendix of 
Ref. 5. We recently obtained the third moment in a six-minute run on 
the 7090. This involved solving a triangular linear system of nine 
equations in nine unknowns. The equations, expressed as polynomials 
in the nine unknowns and five additional parameters, have over 900 
terms. The result is a rational function of the five parameters with a 
numerator of 200 terms and a denominator of 39 terms. The coefficients 
of largest magnitude are 1896 in the numerator and 1460 in the de- 
nominator. The degrees are 1, 1, 3, 7, and 9 for the numerator and 
0, 0, 0, 7, and 9 for the denominator. 

IV. INFORMATION FOR THE ALPAK PROGRAMMER 

4.1 Introduction 

This section is an extension of Section III of Ref. 1. The use and 
implementation of the ALPAK polynomial operations are described 
there, while the use and implementation of the rational-function opera- 
tions are described here. The loading instructions are unchanged except 
that an additional binary deck, called ALPAK3, must be included for 
a run. 

A nonpolynomial rational function is stored as an ordered pair of 
polynomials, namely its numerator and denominator, as illustrated in 
Fig. 1. It consists of a pointer, a rational-function heading, and two 
polynomials stored in the usual way (see Fig. 2 of Ref. 1). The rational- 
function heading contains pointers to the polynomials, which must 
have a common format. 

Integers and polynomials are always recognized as special cases of 
rational functions. If a rational-function operation is used where a 
polynomial operation might have been used, the only penalty will be a 
fraction of a millisecond of additional overhead. 

A rational function can be constructed from its numerator and de- 
nominator polynomials by using RFNDIV (divide) or RFNFRM 
(form). RFNDIV duplicates the two polynomials and constructs the 
rational function from the copies. RFNFRM constructs the rational 
function from the given polynomials and clears their pointers. RFNFRM 
has an optional argument which can be used to indicate that the numera- 
tor and denominator are known to be relatively prime. 

4.2 Input-Output Operations 

RFNRDF F read format (a) 

F RFNCVF (X,15,Y,21,Z,36) convert format (b) 



ALPAK. SYSTEM 



795 



PROGRAM STORAGE 


DATA BUFFER 






R 


POINTER 






RATIONAL- 
FUNCTION 
HEADING 














































NUMERATOR 
HEADING 




DENOMINATOR 
HEADING 










1 


























F 


FORMAT 


! 


NUMERATOR 
DATA 






DENOMINATOR 
DATA 






t 





















Fig. 1 — A rational function R with format F. 



RFNRDD 


R,F 


read data 


(c)* 


RFNCVD 


R,F,HNJID 


convert data 


(d) 


RFNCLR 


R 


clear 


(e) 


RFNSTZ 


R 


store zero 


(0 


RFNSTI 


R 


store identity 


(g) 


RFNSTC 


R,A,B 


store constant 


(h) 


RFNSTV 


R.X.F 


store variable 


(i) 


RFNPRT 


R,CC,(NAME) 


print 


(J) 


RFNPCH 


R,(NAME) 


punch 


(k) 1 


RFNPRP 


R,CC,(NAME) 


print and punch 


(1)' 


RFNRDP 


R,F,CC,(NAME) 


read and print 


(m) 


RFNCVP 


R,F,HN,HD,CC,(NAME) 


convert and print (n) 


A = numerator of constant 






R = denominator of constant; if omittec 


, the denominator 


is un- 


derstood to be one 







* RFNRDD reads two polynomials from cards, interpreting the first as the 
numerator and the second as the denominator of a rational function. If a poly- 
nomial is to be read by RFNRDD, a unit denominator with a complete set of 
zero exponents must be provided. A rational function with a constant numerator 
(denominator) punched by RFNPCH or RFNPRP cannot be read by RFNRDD 
unless a complete set of zero exponents is added to the numerator (denominator). 
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CC = control character for printer 

F = format (symbolic address of /for format statement) 
HN = Hollerith data for numerator (symbolic address of data) 
HD = Hollerith data for denominator (symbolic address of data); 
if omitted, the denominator is understood to be one 
NAME = alternative name for rational function (not exceeding 21 
characters) 
R = rational function (symbolic address of pointer) 
X = variable (specified in the manner indicated by the last previ- 
ous VARTYP declaration. See Ref. 1, Section 3.5). 

(a) RFNRDF F 
Same as POLRDF. See Ref. 1, Section 3.2. 

(b) F RFNCVF (X,15,Y,21,Z,36) 
Same as POLCVF. See Ref. 1, Section 3.2. 

(c) RFNRDD R,F 

Read the rational function R from cards according to format F. R is the 
address of a pointer for the rational function and F is the address of a 
format statement. R must consist of a polynomial numerator and poly- 
nomial denominator punched in cards in that order as specified by Ref. ] , 
Section 3.2. 

(d) RFNCVD R,F,HN,HD 

Same as RFNRDD except that the numerator and denominator poly- 
nomials are to be found in core in blocks of no more than 12 BCI words 
each starting at HN and HD, respectively. 

(e) RFNCLR R 

Clear the rational function R. This clears both numerator and denomi- 
nator polynomials as well as the R heading if R is not itself a polynomial. 
If the R pointer contains zero or points to an idle heading, then RFNCLR 
is a no-op. 

(f) RFNSTZ R 
Same as POLSTZ. See Ref. 1, Section 3.2. 

(g) RFNSTI R 
Same as POLSTI. See Ref. 1, Section 3.2. 
(h) RFNSTC R,CN,CD 
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Same as POLSTC if CD is omitted. See Ref. 1, Section 3.2. If CD is 
present, then CN and CD are the addresses of constants which become 
the numerator and denominator, respectively, of the rational constant 
R = CN/CD. 

(i) RFNSTV R,X,F 

Same as POLSTV. See Ref. 1, Section 3.2. 

(j) RFNPRT R,CC,(NAME) 

Print the rational function It using CC for the control character of the 
first line of print and NAME (not more than 21 characters) for the name. 
If NAME is not provided, "R" will be used for the name; and, if CC is not 
provided, a minus (for triple spacing) is used for the control character. A 
rational function not a polynomial is printed by printing the name on the 
first line, followed by two polynomial prints with the names "NUMER- 
ATOR" and "DENOMINATOR" respectively 

(k) RFNPCH R,(NAME) 

Punch the rational function R on cards using NAME (no more than 21 
characters) for the name. If NAME is not provided "R" will be used for 
the name. A rational function not a polynomial will be punched as two 
polynomials, numerator and denominator in that order. 

(1) RFNPRP R,CC,(NAME) 

Same as RFNPRT followed by RFNPCH. 

(m) RFNRDP R,F,CC,(NAME) 

Same as RFNRDD followed by RFNPRT. 

(n) RFNCVP R,F,HN,HD,CC,(NAME) 

Same as RFNCVD followed by RFNPRT. 

4.3 Arithmetic Operations 



RFNADD 


R,P,Q 


R 


= p + Q 


add 


(a) 


RFNSUB 


R,P,Q 


h' 


= P -Q 


subtract 


(b) 


RFNMPY 


R,P,Q 


H 


= P*Q 


multiply 


(c) 


RFNDIV 


R,P,Q 


R 


= P/Q 


divide 


(d)* 


RFNSST 


G,F(LISTR> 
(LISTV) 


a 


= F( LISTV 
= LISTR) 


substitute 


(e) 



* RFNDIV can form the quotient of any two rational functions provided the 
divisor is not zero. In contrast, POLDIV has a "no divide" return which is used 
whenever the quotient is not a polynomial. 
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Q = dP/dX 


differentiate 


(f) 


Skip iff R = 


zero test 


(g) 


Skip iff R * 


nonzero test 


(h) 


Skip iffP = Q 


equality test 


(i) 


Q = P 


duplicate 


(J) 


R= -R 


change sign 


(k) 



RFNDIF Q,P,X 

RFNZET R 

RFNNZT R 

RFNEQT P,Q 

RFNDUP Q,P 

RFNCHS R 

F,G,P,Q,R = rational functions (symbolic addresses of pointers) 

X = variable (specified in the manner indicated by the last 
previous VARTYP declaration) t 
LISTR = list of rational functions. 

LISTV = list of variables (specified in the manner indicated by the 
last previous VARTYP declaration)! 



4.3.1 Nutation 

In the following descriptions, if R denotes a rational function, we de- 
note its numerator by RN and its denominator by RD. In particular, if 
R is a polynomial, RD is the constant polynomial 1. 

4.3.2 Descriptions 

(a) RFNADD R,P,Q 

The inputs P and Q are rational functions in canonical form. First 
POLGCD is used to obtain 

G = (PD,QD). 

Then the polynomials 

AN = PN*(QD/G) + QN*(PD/G) 
and 

AD = (PD/G)*(QD/G)*G 

are computed. Note that the parenthesized fractions are polynomials. 
Next, 



is obtained and 



H = (AN,G) 
(AN/H) 



R = 



(AD/H) 



t See Ref. 1, Section 3.5. 
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is formed. RN and RD are now relatively prime (the proof of this is an 
exercise for the reader), so POLGCD need not be used when placing R 
in canonical form. 

(b) RFNSUB R,P,Q 

RFNCHS [see (k) below] and RFNADD are applied to compute 

R = P + (-Q). 

(c) RFNMPY R,P,Q 
First the functions 

^1 = PN/QD 

and 

B = QN/PD 

are formed and placed in canonical form. Then 
R = AN*BN/AD*BD 

is formed. RN and RD are already relatively prime, so POLGCD need 
not be used when placing R in canonical form. 

(d) RFNDIV R,P,Q 

RFNDIV is identical to RFNMPY except that Q must be not zero and 
the roles of QN and QD are interchanged. If Q is zero, the diagnostic re- 
mark "ZERO DENOMINATOR" is printed and the job is termi- 
nated. 

(e) RFNSST G,F(LISTR) (LISTV) 

RFNSST is exactly the rational function equivalent of POLSST; in par- 
ticular, the format constraints on F and G are identical. If F is a poly- 
nomial, the rational functions of LISTR are substituted for the varia- 
bles in LISTV term-by-term to accumulate the final result. If F is a 
rational function, this procedure is applied to the numerator and de- 
nominator polynomials of F in succession and the resulting rational func- 
tions are divided (using RFNDIV) to obtain G. 

(f) RFND1F Q,P,X 
First POLDIF is used to compute 

PN' = d(PN)/dX 
PD' = d(PD)/dX, 
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and POLGCD is used to obtain 

G = (PD,PD'). 

Next the polynomials 

AN = (PD/G)*PN' - PN*(PD'/G) 
and 

AD = PD*(PD/G) = (PD/G) 2 *G 

are computed. Note that the parenthesized fractions are polynomials. 
Finally 

H = (AN,G) 
is obtained, and 

(AN/H) 



Q = 



(AD/H) 



is formed. Since QN and QD are relatively prime, POLGCD need not 
be used when placing Q in canonical form. 

(g) RFNZET R 

Same as POLZET. See Ref. 1, Section 3.3 

(h) RFNNZT R 

Same as POLNZT. See Ref. 1, Section 3.3. 

(i) RFNEQT P,Q 

If P and Q are both polynomials, POLEQT is applied. If only one of 
them is a polynomial, they are not equal. If neither is a polynomial, 
POLEQT is applied to both numerators and both denominators. If the 
rational functions are unequal, the next instruction is executed, if they 
are equal, then the next instruction is skipped. 

(j) RFNDUP Q,P 

Q is replaced by a copy of P. 

(k) RFNCHS R 

If R is a polynomial, POLCHS is applied to R. If R is not a polynomial, 
POLCHS is applied to its numerator polynomial . 

4.4 Truncated Power Series Operations 

ALPAK contains two macros for dealing with truncated power series 
with rational function coefficients. These are RFNTRC (truncate) and 
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RFNMPT (multiply and truncate). Addition can be handled with 
RFNTRC and RFNADD. Each truncated power series must be stored 
as a rational function in a format whose first k variables are the power 
series variables. The denominator, if any, must be independent of these 
variables. The command 

RFNTRC P,ORD,K 

where K contains the number of power series variables k and ORD con- 
tains an integer n, causes P to be truncated to order n. That is, all terms 
of order greater than n in the first k variables are deleted. The command 



R,ORDR,P,ORDP,Q,ORDQ,K 



RFNMPT 

is represented by the equation 

R = P*Q 

where P and Q are truncated power series. K is the address of the num- 
ber of power series variables, ORDP and ORDQ are the addresses of the 
orders of P and Q respectively, and ORDR is an address for the order 
of R, which is to be comupted. If P(Q) contains any terms of order 
greater than ORDP(ORDQ), they will be deleted. 

4.5 Miscellaneous Operations 

(Caution: read descriptions carefully.) 
G,A,B 



POLGCD G,A,B greatest common divisor of (a) 

polynomials 

INTGCD greatest common divisor of in- (b) 

tegers in AC and MQ 

PWVSTO XK,K,W,FA store a power of the TFth vari- (c) 

able 

VARNUM W,X,FA variable number (d) 

RFNFRM R,N,D form (e) 

EXPAND N,D,R,IORP expand (f) 

SUBLCK PJ,P,J,W sub-block (g) 

PWVFAC K,P,W factor off a power of the TFth (h) 

variable 

DEGREE K,P,W degree. (i) 

(a) POLGCD G,A,B 

Replace G by a greatest common divisor of the polynomials A and B. 

(b) INTGCD 
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Replace the integer in the AC by the greatest common divisor of it and 
the integer in the MQ. 

(c) PWVSTO XK,K,W,FA 

Replace XK by the Kth power of the Wth variable in the format whose 
address is at FA. K is the address of the power, and W is the address of 
the variable number. 

(d) VARNUM W,X,FA 

Replace the contents of W by the variable number of the variable X in 
the format whose address is at FA. X is the address of the variable name 
in BCI. 

(e) RFNFRM R,N,D 

Same as RFNDIV except: (i) N and D must be polynomials; (ii) N 
and D become the property of R, and their pointers in the calling pro- 
gram are replaced by zeros; (in) R must be distinct from N and D; 
and (iv), if N and D are known to be relatively prime, a fourth argu- 
ment NOGCD can be added to the calling sequence in order to save the 
time which would otherwise be spent in finding their greatest common 
divisor. 

(f) EXPAND N,D,R,IORP 

This is the inverse of RFNFRM. N and D must initially contain zeros. 
They are filled in with pointers to the numerator and denominator of R 
respectively. The R heading is marked as idle. If R is an integer or a poly- 
nomial, N is filled in with a pointer to R, the R pointer is replaced by 
zero, and control is transferred to IORP. If IORP is omitted, control is 
transferred to the next instruction. 

(g) SUBLCK PJ,P,J,W 

P must be a polynomial independent of the first W — 1 variables, if 
any. Then, by the definition of the polynomial canonical form, the terms 
of P are ordered according to the powers of the Wth variable. SUBLCK 
replaces PJ by the polynomial consisting of that sub-block of P, if any, 
whose terms all involve the Jth power (0 < J < degree of P) of the 
Wth variable. If P contains no terms involving the Jth power of the 
Wth variable, SUBLCK replaces PJ by the zero polynomial. 

(h) PWVFAC K,P,W 

P must be a polynomial independent of the first W — 1 variables, if 
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any. PWVFAC replaces the contents of K by the smallest exponent of 
the Wth variable in P, and divides P by that power of the TFth variable. 

(i) DEGREE K,P,W 

P must be a polynomial independent of the first W — 1 variables, if 
any. DEGREE replaces the contents of K by the degree of P in the 
IPtli variable. 



v. OUTLOOK 

A new version of ALPAK (to be called ALPAKB) is now being 
developed. Its foundation is a programming system (see Ref. 4) called 
STGPAK (storage package), which provides (i) dynamic storage 
allocation, (ii) automatic recursion, and (Hi) "delayed -decision diag- 
nostics." 

The storage allocation orders make it possible to obtain contiguous 
blocks of storage of arbitrary length as needed (provided that sufficient 
space is available) and to return idle space to the system. A block may 
contain sub-blocks and /or pointers to other blocks. This will permit 
the introduction of higher level data structures including formal prod- 
ucts of polynomials, thereby helping to alleviate the greatest-common- 
divisor problem. 

The use of a public push-down list for subroutine storage makes 
recursive programming fully automatic. That is, a subroutine can call 
itself without taking special measures to preserve its arguments and 
intermediate results. The diagnostic facilities permit the decision re- 
garding what to do about an overflow (shortage of space or time) or 
error detected in a given subroutine, to be delayed until control has been 
returned to some higher level subroutine or to the main program. 

The authors hope that STGPAK together wHh a macro compiler 
now being developed by Miss D. C. Leagus and W. S. Brown will 
simplify and expedite the programming of ALPAKB subroutines. The 
compiler should also be useful in the writing of main programs. 

Apart from these matters, which are not directly related to algebra, 
our plans for ALPAKB include multiple precision integer arithmetic, 
an improved strategy for finding greatest common divisors, and a 
complete set of operations for truncated power series. 
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